Jamming versus Glass Transitions 
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Recent ideas based on the properties of assemblies of frictionless particles in mechanical equilibrium provide 
a perspective of amorphous systems different from that offered by the traditional approach originating in liquid 
theory. The relation, if any, between these two points of view, and the relevance of the former to the glass 
phase, has been difficult to ascertain. In this paper we introduce a model for which both theories apply strictly: 
it exhibits on the one hand an ideal glass transition and on the other 'jamming' features (fragility, soft modes) 
virtually identical to that of real systems. This allows us to disentangle the two physical phenomena. 
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The traditional way to introduce the glass transition is to 
start with a liquid of, say, hard particles at low pressures Jl|, 
and to consider a slow compression. At a given point, the 
viscosity increases dramatically, the dynamics becomes slug- 
gish and the system falls out of equilibrium. Slower com- 
pression protocols push the equilibrium regime further and 
make the transition sharper, and one conceives of a limit of 
infinitely slow annealing in which one recognizes (perhaps) 
a true thermodynamic change of state. An argument in favor 
of a thermodynamic transition was given years ago by Kauz- 
mann, who interpreted it as a consequence of the liquid run- 
ning out of configurational entropy. Although the possibility 
of proving the occurrence of such an ideal glass transition for 
a real system seems remote, there is a family of models (or 
approximations), for which this picture holds strictly, within 
the so-called 'Random First Order' (RFO) scenario yiLa], a 
Mean-Field theory allowing for a complete analytic analysis. 

An apparently unrelated set of ideas comes from consid- 
ering amorphous assemblies of hard, frictionless particles in 
mechanical equilibrium. In general, such systems can be hy- 
postatic, hyperstatic or isostatic, depending on whether the set 
of contacts yields a number of conditions smaller, larger, or 
precisely equal to the number of degrees of freedom - just 
like a table with two, four or three legs touching a floor. A 
polydisperse system of spheres will with probability one be 
isostatic much in the same way that a table with 

many legs will only have three touching on a rough floor. 
Because breaking one single contact already destabilizes an 
isostatic system, one can then argue that such systems are 
marginal ItZl, [fa, |9|] : they have large responses and a spectrum 
of vibrations with finite density of very low frequency 'soft' 
modes. One is then in the presence of a "Jammed" configu- 
ration that is critical, with diverging lengths, nontrivial expo- 
nents, etc BSD. 

A natural question is whether this critical "Jamming" phe- 
nomenon is in some way a finite-dimensional (beyond mean- 
field) manifestation of the glass transition |9D or in other 
words, whether the diverging length associated with the soft 
modes of an amorphous packing is a manifestation of the or- 
der underlying the glass transition in finite dimensions. In this 
paper we propose to answer such questions by constructing 
a family of models having both theories written side-by-side 




FIG. 1: The model: Particles evolve in a d-dimensional space (right), 
but are only able to 'see' a pre-established subset of the rest. Which 
particle interacts with which is encoded in the quenched graph (left). 

with a 'J-point' isostatic equilibrium with a spectrum of soft 
modes virtually identical to that of systems of hard spheres 
and that are by construction mean-field models with a 'Ran- 
dom First Order' behavior J^,[3|]. 

The model — Consider the usual hard-sphere model 
where N non overlapping spheres live in a cube of dimension 
d, with periodic boundary conditions. Following an approach 
introduced in 1 1 Oil - we now obtain a mean field version of this 
model: the idea is to let each sphere interact only with a set 
of z other ones, all the rest being transparent to it. The set of 
particles that interact with a given sphere is chosen once and 
for all. If z = N — 1 for each particles, we recover the original 
problem, but we shall instead work with z finite. The model is 
completely defined by a quenched regular random graph (see 
FIG[T|) with nodes labeled by the particle number and bonds 
that denote interaction. 

This model presents a number of advantages: first, since 
interactions are set through a random tree-like graph — or 
Bethe-Lattice — it is by construction a mean-field glass model 
that can be studied exactly with the cavity method [11] which 
is an extension of the well-known work of Bethe and Peierls 
on trees. Secondly, our model can be seen as a constraint satis- 
faction problem (CSP) defined on a random graph, where the 
constraint on each node is that the position within the cube 
of the sphere on that node be such that it does not overlap 
with any of the other z interacting spheres (that are linked on 
the graph). This model bridges the gap between the field of 
random discrete CSP such as the coloring problem 111 2i 1 1311 
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and disordered hard-sphere packings ifioll . Our approach thus 
provides an original way to attack packing problems within 
mean field theory. In this paper, we take the first steps in this 
direction by studying numerically the model in d = 2 and re- 
strict ourselves to the question of determining the differences 
between the Jamming and the glass transition. 

Glassy behavior — The low pressure/low density "liq- 
uid" state can be discussed easily within the cavity method. 
It just corresponds to considering the model on a tree in the 
usual Bethe-Peierls way while ignoring long range correla- 
tions. A straightforward computation shows that the entropy 
of the liquid state on a random graph of connectivity z is 
Suquid = z log (1 — ttD 2 )/2, where D is the diameter of the 
spheres IH4J1. Using a pressure — P conjugated to the sphere 
volume llOll V = D 2 this yields the equation of state in the 
liquid phase P = (z7r/2)(l — irD 2 )^ 1 . However, this equa- 
tion is inconsistent at large pressures: indeed for P = oo 
it yields a limiting value independent of the connectivity z, 
which is unphysical, and where each sphere is in contact with 
all its neighbors; this requires not only the number of con- 
tacts to be much larger than the isostatic value, but also the 
presence of a 2— coloring of the graph 1 10] and thus violates 
rigorous results in graph theory: this proves that a phase tran- 
sition must intervene at a finite pressure value. It is indeed 
possible to prove within the cavity method that the liquid 
phase is unstable towards a glass phase at some finite pres- 
sure and ongoing analytic work in this direction shows that 
the system has a glass phase lfl4ll of the Random First Or- 
der kind for z sufficiently large, as usual with frustrated mod- 
els on such structures till 
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is illustrated on Figj2] in a system with a slight polydisper- 
sity 01811 (in which case one should compute the averaged 
Suquid — z\og (1 — 7r(i?i + i?2) 2 )/2 over the distribution of 
radii P(R)): a clear signal of glass transition with its charac- 
teristic compression rate dependence is observed. 

Jamming — We now turn to jammed configurations ob- 
tained following the procedure of |9[] : starting from small par- 
ticles we 'inflate' them infinitesimally, and adjust their posi- 
tions to eliminate any (infinitesimal) overlap, until a jammed 
state is obtained. At very large pressure, particles are found to 
be either locally blocked by their neighbors (in which case 
they have at least three contacts), or are 'rattlers' that can 
be displaced without moving the others 0. The locally 
blocked particles have a subset that cannot be displaced even 
with collective rearrangements - except global translations. 
For our polydisperse system, this subset turns out not to have 
redundant contacts, and hence constitutes an isostatic core 
(see left panel of FIG. |3}. We find many states with slightly 
different densities, all of which are isostatic and have indis- 
cernible properties. 

The jammed configurations have a spectrum of modes that 
is strikingly similar to those found in the ordinary finite di- 
mensional system Ji, 2(1 21 ]. In FIG.|3]we show the vibration 
modes obtained by considering a soft potential and a small 
superposition of particles that follows from inflating the par- 
ticles slightly beyond the jammed configuration, as in Ref. 
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FIG. 2: Volume versus inverse pressure for a system of N = 4000 
particles with z = 100 and 9% polydispersity. The choice of vari- 
ables for the axes is to stress the analogy with the more usual energy- 
temperature annealing plots. The annealing curves follow the liquid 
analytical solution at low pressures, but break away at the glass tran- 
sition pressure P g , which is clearly visible. 



H20tl . A second, perhaps more satisfactory way to study the 
vibrations of a hard sphere ensemble is to consider the very 
high pressure dynamics of hard particles around a jammed 
state, and computing the displacement correlations (Aij) with 
Aij = Xi(t)xj(t) — (xi)(xj) where (•} denotes average over 
a time that is long yet insufficient for escaping the vicinity 
of the jammed configuration. The displacements fall into two 
classes: the many-particle vibrations that scale as P _1 , and 
the rattlers that do not, because the cages are roughly indepen- 
dent of the pressure. The role of the frequencies in a potential 
is played by the eigenvalues of P^ 1 A^ 1 / 2 . The spectrum of 
these eigenvalues is shown in FIG. [3] at infinite pressure there 
is a proliferation of soft modes, but as the volume fraction <p is 
lowered, and the pressure becomes finite, the system interacts, 
due to collisions and rattling, with more particles than what 
is imposed by the isostaticity condition, and the isostaticity- 
related soft modes start disappearing. Again, this is in perfect 
agreement with the usual particle systems [SJ]. 

Jamming versus Glass — Having discussed the isostatic 
packings with the standard J-point phenomenology lIsL foL [2oL 
2211 . as well as the mean field glassy nature of our model, we 
can now distinguish the 'Jammed' and the 'Glass' features as 
follows (FIG.IU). There is a liquid equilibrium line, which ter- 
minates at the equilibrium glass transition pressure Pk- Be- 
low Pk the Gibbs measure is dominated by a few of the deep- 
est glassy states [4]. There are also metastable ones (indicated 
by the horizontal lines), which are stable (only) within the 
mean-field picture. As in any system with a dynamic tran- 
sition, any compression process ending at infinite pressure 
leads the system into a metastable state. The 'jamming line' 
{P = oo , ip < (p < (pj} is the set of such (out of equi- 
librium) blocked configurations. The slower the compression, 
the denser the target state is, and one needs an infinitely slow 
process to reach ipo. In particular, we denote the J point (pj 
as the result of the fastest compression starting from a random 
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FIG. 3: Left: Integrated pair correlation function restricted to the collectively jammed core. The average number of neighbors is four, as one 
expects from an isostatic system in two dimensions. The inset shows a zoom in the small distance region, the horizontal line is the isostaticity 
condition, while the vertical line is the gradient descent step used to jam the system. Center: Vibration modes with small overlaps (here noted 
D — Dj, difference between actual diameter of the particles and diameter particles should have to be at jamming). The resulting spectrum is 
almost identical to that of a particle system (cf FIG. 1 of ref. [ 20]). Right: Hard particles: vibration modes of the isostatic core, rescaled with 
the pressure. In the inset, gap in the spectrum vs. pressure scales as P _1,/2 . The curves are virtually identical to those of a finite dimensional 
system, as obtained by Brito and Wyart 1 8]. These curves correspond to a configuration with N=1000 and z=13. 



configuration (9|], which we expect to be different from the 
deepest level (p a - just as a quench to zero temperature does 
not terminate in the ground state of any complex system. In a 
finite dimensional system, only two lines in the phase diagram 
are stable, apart from the liquid state: the jamming line (at in- 
finite pressure) and the equilibrium glass phase - if it exists at 
all. Their only common point is ((p Q , P = oo). 

A measure of the criticality of a state is the staggered dis- 
placement (£\ ) I h produced by random forces /j = h^. L , 
with £j random unit vectors. This corresponds, in the spin- 
glass literature, to the quantity (1 — qea)/T. It is proportional 

to xea = P 2 £,«*?(*)_> - (*i) 2 ) = P 2 M(A)). If lo are 
the eigenvalues of P~ 1 A~ 1 / 2 and Q{uj) the density of modes, 
then xea ~ / dujQ(ui)/uj 2 . Due to the soft modes BSD, this 
quantity diverges on approaching the jamming line when nor- 
malised this way, see inset of FIG. [3] A perhaps more standard 
measure is the function Xi(t, t') = • Aij(t)Aij(f). It has 
a maximum which diverges as the system approaches either 
the jamming points (because of the soft modes) or the equilib- 
rium glass line, because of the divergence of activation time 
(i.e. the time needed to overcome a free-energy barrier). 

Thinking in terms of landscape iflOD . the picture that 
emerges is one of a system with many glassy states, separated 
by high and wide free-energy barriers, as usual in RFO mod- 
els. Fragility and soft modes are then a property of the (infinite 
pressure) jammed configurations within each of these states, 
but states are not marginal in the coarse grained view obtained 
as soon as temperature is nonzero and the pressure finite (at 
least in all but the more superficial 'threshold' levels). This 
separation between large valleys with critical bottom is con- 
sistent with the fact that assemblies of frictionless particles at 
zero temperature are at the same time fragile (even slight shear 
stresses make them creep some amount) while still capable of 




FIG. 4: A sketch of the density(<^)-pressure(P) plane in RFO models 
(see e.g. |3, 10, 16]). Dynamic and static transition pressures are Pd 
and Pk- The glass phase and the 'jamming line' (P — oo,ip < 
ip < ifij) are clearly distinct. The isostaticity-related quantity Xea 
is finite within a state if P is finite (except perhaps for a small effect 
of acoustic modes in d — 2). In particular it is finite within the ideal 
glass state, while it is infinite on the jamming line. Everywhere in the 
trapezoidal region delimited by the equilibrium line < P < Pd, 
the threshold level, and the jamming line the activation time diverges 
(at times exponential in N), and with it the four-point function Xi- 
In finite dimensions Xi would only diverge strictly at the jamming 
line (at finite times) and close to the glass line (at times comparable 
to the activation time). 



resisting a shear stress proportional to the pressure without 
flowing continuously (they have internal friction [23].) 

Jamming a Crystal — To unambiguously demonstrate that 
a large stable state with an internal 'fragile' structure is pos- 
sible in finite dimensions, and that the Jamming point is of 
different nature than that of the glass transition, let us finally 
consider an ordinary FCC crystal of spherical particles with 
very small (.003%) polydispersity. The crystalline order is 
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preaching the glass line from the 
ing the jamming line at P — oo B24tl 



liquid phase, than approach- 



FIG. 5: Left: A jammed configuration in the polydisperse crystal 
(N=864). The distortion of the crystalline order due to polydispersity 
is extremely weak. Right: Spectrum of normal modes of the poly- 
disperse crystal, obtained in the same way as FIG.3, central panel. 
There are many soft modes (To be compared with FIG. [3}. In inset: 
integrated pair correlation function, including and excluding rattlers: 
the number of contacts of non-rattlers is m 6. 



hardly affected (see FIG. [3), and polydispersity becomes ir- 
relevant at finite pressures P. However, at P = oo the system 
is isostatic, and, not surprisingly, has a spectrum of soft modes 
even richer than that of an amorphous packing. In an amor- 
phous solid, each one of the many equilibrium glass states 
plays a similar role to the one of the crystal state above, their 
very high pressure jamming properties being juxtaposed with 
the underlying long-range glass order, already established at 
finite pressure. 

Conclusion — In this paper we confront two different vi- 
sions of amorphous systems: a glassy solid state with order 
characterized by a permanent, amorphous modulation of den- 
sity - not unlike a crystal or a quasi-crystal - and a jamming 
situation brought about by chains of force associated with ac- 
tual contact between hard particles. These two phenomena 
may coexist but are distinct. The fact that a mean-field model 
reproduces both the glass transition and J-point criticality in a 
separate way suggests that we abandon the idea that the for- 
mer is some kind of finite-dimensional realization of the latter. 
The mean-field picture just above the glass transition pressure 
is one of large basins separated by high barriers, without an 
excess of truly zero frequency modes except perhaps at the 
least deep states. The isostaticity-related marginality of the 
jammed configurations appears only at very large pressures, 
deep within a basin. This local criticality combined with ab- 
sence of criticality 'in the large' is attested by the paradoxical 
fact that amorphous matter is fragile to small stresses but may 
still sustain extensive stresses without flowing, a fact that can 
be understood easily in the example of the polydisperse crys- 
tal. Perhaps the most evident manifestation of the different 
nature of the jamming (P = oo) and the glass lines is the fact 
that, although lengths such as associated to four-point function 
\4 diverge in both, the growth is astronomically slower ap- 



In conclusion, we have introduced a set of new models that 
can be studied analytically and numerically and that provide 
non-trivial connections between different fields. We hope this 
will be useful to attack the glass and granular problems ana- 
lytically in the future. 
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with us their unpublished results 11411 . 
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